In [34]:
%load_ext autoreload
%autoreload 2


The autoreload extension is already loaded. To reload it, use:
  %reload_ext autoreload

In [35]:
import matplotlib.pyplot as plt
%matplotlib inline
from plotter import StackHist, StackHistWithSignificance
from utils import ResultSet

In [36]:
rs_TTZ  = ResultSet("TTZ",  "../data_ichep/TTZToLLNuNu_treeProducerSusyMultilepton_tree.root")
rs_TTW  = ResultSet("TTW",  "../data_ichep/TTWToLNu_treeProducerSusyMultilepton_tree.root")
rs_TTH = ResultSet("TTH", "../data_ichep/TTHnobb_mWCutfix_ext1_treeProducerSusyMultilepton_tree.root")
rs_TTTT = ResultSet("TTTT", "../data_ichep/TTTT_ext_treeProducerSusyMultilepton_tree.root")


Loading unchanged result file  ../data_ichep/TTZToLLNuNu_treeProducerSusyMultilepton_tree_result.root
Loading unchanged result file  ../data_ichep/TTWToLNu_treeProducerSusyMultilepton_tree_result.root
Loading unchanged result file  ../data_ichep/TTHnobb_mWCutfix_ext1_treeProducerSusyMultilepton_tree_result.root
Loading unchanged result file  ../data_ichep/TTTT_ext_treeProducerSusyMultilepton_tree_result.root

In [37]:
def stack(attr, ax, title, cls=StackHist, stack_signal=False, scale = None, **kwargs):
    hs = cls(title)
    hs.add_mc_background(getattr(rs_TTZ, attr), 'TTZ', lumi=rs_TTZ.lumi, plot_color='r')
    hs.add_mc_background(getattr(rs_TTW, attr), 'TTW', lumi=rs_TTW.lumi, plot_color='g')
    hs.add_mc_background(getattr(rs_TTH, attr), 'TTH', lumi=rs_TTH.lumi, plot_color='b')
    if scale is None:
        scale = 1 if stack_signal else 20
    hs.set_mc_signal(getattr(rs_TTTT,attr), 'TTTT', lumi=rs_TTTT.lumi, stack=stack_signal, scale=scale, plot_color='k')
    
    hs.luminosity = 40
    hs.energy = 13
    hs.logy=False
    hs.draw(ax, **kwargs, save=True)

The main discriminating variables in the analysis are Jet and B-Jet Multiplicities

Below are shown the distributions for the TTTT as well as relevant BG processes. The "baseline" selection for the trilepton channel is:

  • Exactly three leptons
  • Three or more B-Jets
  • Z-Mass Window veto: The invariant mass of all same-flavor opposite-sign lepton pairs must be outside the Z-Mass window (70,105)GeV

In [38]:
fig, ((ax0, ax1), (ax2, ax3)) = plt.subplots(2,2, figsize=(26,20), tight_layout=False)
stack('jet_count', ax0, 'Jet Multiplicity', cls=StackHistWithSignificance, high_cut_significance=True)
stack('jet_count_trilepton', ax1, 'Jet Multiplicity ($N_{lep}=3$)', cls=StackHistWithSignificance, high_cut_significance=True)
stack('jet_count_b_jet3', ax2, 'Jet Multiplicity ($N_{Jet_B}\ge3$)', cls=StackHistWithSignificance, scale=1, high_cut_significance=True)
stack('jet_count_z_mass_veto', ax3, 'Jet Multiplicity (Z-mass Veto)', cls=StackHistWithSignificance, scale=1, high_cut_significance=True)



In [39]:
fig, ((ax0, ax1), (ax2, ax3)) = plt.subplots(2,2, figsize=(26,20), tight_layout=False)
stack('b_jet_count', ax0, 'B-Jet Multiplicity', cls=StackHistWithSignificance, high_cut_significance=True)
stack('b_jet_count_trilepton', ax1, 'B-Jet Multiplicity ($N_{lep}=3$)', cls=StackHistWithSignificance, high_cut_significance=True)
stack('b_jet_count_b_jet3', ax2, 'B-Jet Multiplicity ($N_{Jet_B}\ge3$)', cls=StackHistWithSignificance, scale=1, high_cut_significance=True)
stack('b_jet_count_z_mass_veto', ax3, 'B-Jet Multiplicity (Z-mass Veto)', cls=StackHistWithSignificance, scale=1, high_cut_significance=True)



In [40]:
fig, ((ax0, ax1), (ax2, ax3)) = plt.subplots(2,2, figsize=(26,20), tight_layout=False)
stack('jet_count', ax0, 'Jet Multiplicity', cls=StackHistWithSignificance, high_cut_significance=True)
stack('jet_count_trilepton', ax1, 'Jet Multiplicity ($N_{lep}=3$)', cls=StackHistWithSignificance, high_cut_significance=True)
stack('jet_count_eem_trilepton', ax2, 'Jet Multiplicity ($e^+e^-\mu$)', cls=StackHistWithSignificance, high_cut_significance=True)
stack('jet_count_mme_trilepton', ax3, 'Jet Multiplicity ($\mu^+\mu^-e$)', cls=StackHistWithSignificance, high_cut_significance=True)



In [41]:
fig, ((ax0, ax1), (ax2, ax3)) = plt.subplots(2,2, figsize=(26,20), tight_layout=False)
stack('b_jet_count', ax0, 'B-Jet Multiplicity', cls=StackHistWithSignificance, high_cut_significance=True)
stack('b_jet_count_trilepton', ax1, 'B-Jet Multiplicity ($N_{lep}=3$)', cls=StackHistWithSignificance, high_cut_significance=True)
stack('b_jet_count_eem_trilepton', ax2, 'B-Jet Multiplicity ($e^+e^-\mu$)', cls=StackHistWithSignificance, high_cut_significance=True)
stack('b_jet_count_mme_trilepton', ax3, 'B-Jet Multiplicity ($\mu^+\mu^-e$)', cls=StackHistWithSignificance, high_cut_significance=True)